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Abstract 

We present here a theoretical model in order to describe the competi- 
tion between the Kondo effect and the spin glass behavior. The spin glass 
part of the starting Hamiltonian contains Ising spins with an intersite ex- 
change interaction given by the local van Hemmen model, while the Kondo 
effect is described as usual by the intrasite exchange Jk- We obtain, for 
large Jk values, a Kondo phase and, for smaller Jk values, a succession, 
with decreasing temperature, of a spin glass phase, a mixed spin glass- 
ferromagnetic one and finally a ferromagnetic phase. This model improves 
the theoretical description of disordered Kondo systems with respect to 
previous models and can account for experimental data in Cerium disor- 
dered systems like CeCui- x Ni x alloys. 

1 Introduction 

The interplay between disorder and strong electronic correlations is recognized 
as a very interesting issue in condensed matter physics. There are now many 
experimental evidences showing the very important role of the disorder in f- 
electron systems in addition to the RKKY or Kondo interactions pQ. As a result, 
it can appear complex phase diagrams which show spin glass (SG) phases in ad- 
dition to the onset of antifcrromagnetism (AF) or ferromagnetism (FE), regions 
dominated by the Kondo effect, the presence of Quantum Phase Transitions 
(QPT) and exotic regions which present non-Fermi liquid behavior (NFL) 2 . 

Earlier experimental results can illustrate the mentioned complexity. For 
instance, in CeAu\- x Co x Siz alloys 13!, when Au is replaced by Co, it first 
appears a SG phase, then there is the onset of an AF phase with the Neel 
temperature decreasing towards a Quantum Critical Point (QCP). Thus, the 
glassy behaviour tends to decrease with the increase of x and finally, for a; > 0.9, 
there is a complete screening of magnetic moments due to the Kondo effect. 

More recently, experimental findings in CePd\- x Rh x [H[S] and CeNi\- x Cu x 
[5J [7] have enlarged the set of non-trivial behaviour in disordered /-electron sys- 
tems. In both systems, there are strong indications that a glassy behaviour is 



present in a suitable range of doping and this behaviour has been recently iden- 
tified to appear as a cluster glass state. In the well studied CeNi\- x Cu x case, 
the Kondo interaction is dominating for x smaller than approximately 0.2 [5]. 
However, the intermediate doping regime has been extensively studied both ex- 
perimentally and theoretically and finally a complex scenario is obtained when 
the temperature is decreased. In the first experimental studies on CeNi\- x Cu x 
alloys with x typically between 0.3 and 0.6, a SG phase has been obtained below 
the paramagnetic state and then there is a transition to a ferromagnetic phase 
at lower temperatures. 

More sophisticated experiments have recently shown that dynamic magnetic 
clusters are developping at low temperatures below the paramagnetic state. 
More precisely, there is the formation of clusters due to short range ferromag- 
netic correlations below a certain temperature T* . The volume fraction of these 
clusters increases as temperature is lowered and they become frozen at T c i well 
below T* and, therefore, it appears an inhomogeneous ferromagnetic order at 
very low temperatures [3 [7]. Thus, there is a change, below the paramagnetic 
phase, from a cluster spin glass to a disordered ferromagnetic order without any 
sharp transition, but with a mixed and disordered intermediate phase. 

A Kondo-Cluster-Glass state has been also recently evidenced in CePd\- x Rh 
alloys at very low temperatures. This system exhibits a continuous evolution 
from a ferromagnetic order in CePd, with a Curie temperature T c — 6.6K , to an 
intermediate-valence ground state in CeRh. The Curie temperature decreases 
continuously with increasing x and tends to 25 mK at the value x = 0.87. De- 
spite pronounced non-Fermi-liquid behavior in the proximity of this concentra- 
tion for specific heat and thermal expansion, it was concluded from the analysis 
of the Gruneisen ratio that there is no QCP [9]. On the opposite, a "Kondo- 
cluster-glass" state was found for x larger than 0.65: there is firstly the formation 
of clusters with predominantly ferromagnetic couplings of the /-moments below 
a given temperature T* and then a random freezing of the cluster moments be- 
low a smaller temperature T c i [H[5]. Thus, there are clearly similarities between 
the low temperature behaviors of CeNi\- x Cu x and CePd\- x Rh x alloys, but 
both a more profound analysis of the different data and the role of the Kondo 
effect have to be precised in these two systems. 

Several theoretical studies have tried, since already some time, to account for 
the previous experimental data. A Kondo lattice with an additional Ising term 
and a random coupling between localized spins, called here the Kondo-Ising 
Lattice (KIL) model [TOj, [TTJ [12] , has been firstly used to study the competi- 
tion between the Kondo effect and magnetism when disorder is present within 
the Static Approximation (SA) [13]. It appears that, for CeAui^ x Co x Si^ al- 
loys, a Gaussian random distributed bond would be adequate as can be seen 
in Refs. [T21 HI]. The same model has been also firstly used to describe the 
case of CeNii_ x Cu x alloys, where the disorder has been introduced within the 
classical Sherrington-Kirkpatrick (SK) model [15] by taking Gaussian random 
intersite coupling Jy with a mean value Jo different from zero to describe the 
ferromagnetic ordering The phase diagram giving the temperature T ver- 
sus the strength J# of the Kondo interaction has been computed and we have 
obtained, besides the Kondo state, magnetic phases like Spin Glass (SG), Fer- 
romagnetic (FE) and a mixed phase (SG+FE). For this particular solution, the 
ferromagnetic order occurs with replica symmetry breaking. This phase dia- 
gram could be, therefore, a good starting point to describe the scenario found 
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in CeNii_ x Cu x alloys. Unfortunately, for this particular kind of disorder, the 
Curie temperature T c is always higher than the freezing one, which is a scenario 
opposite to the experimental situation observed in CeNi\- x Cu x system. 

Thus, in order to solve the preceding difficulty, a completely different per- 
spective has been adopted in reference [16]. The theoretical description of the 
disorder has been modified from a bond disordered coupling to a site disordered 
one. In that case, the Jij coupling is a generalization of the Mattis model [T7] 
used extensively to study complex systems [TS], given as Jij = ^ Su=i £f £j> 
where £f is a random variable which follows a bimodal distribution. 

One important aspect is that, in the corresponding mean field approach using 
such values, it is possible to introduce a parameter which allows to control 
the level of frustration in the problem |16j . The first interesting result is that 
the Kondo solution is robust in the large Jk limit, no matter what is the level 
of frustration. For weak frustration and small Jk , below a certain temperature, 
it appears a SG solution. When the temperature is further decreased, the SG 
solution is replaced by Mattis states which have the same thermodynamics as 
a FE phase [TS]. This result suggests that the situation found in CeNi\_ x Cu x 
alloys would be an example of weak frustration. Nevertheless, there is an impor- 
tant difference between this model and the approach of the Gaussian distributed 
Jij. For the kind of disorder given by this generalized Mattis model, there is 
no mixed phase solution for the order parameters, but on the contrary, there is 
a first order phase transition between the SG and Mattis states; such solutions 
can obviously coexist, but one of them is always metastable. In conclusion, 
our previous Mattis-like model gives the SG phase above the FE phase, but it 
cannot yield a real SG+FE mixed Phase |16| . 

Thus, in order to improve the preceding description and to have, therefore, 
a better agreement with experiment, we introduce here, in our previously used 
KIL model, a new kind of site disordered coupling Jy, originally introduced 
by van Hemmen (vH) to study the Spin Glass in the classical Ising model [TP] , 
The phase diagram obtained from such a classical model displays not only SG, 
FE+SG and FE phases, but also they can appear in that order when temper- 
ature is decreased. In this particular case, the SG+FE phase is characterized 
by both non zero magnetization and SG order parameters. Recently, a work 
[2D] has studied a mean field solution of a quantum version of the vH model 
with an applied transverse field T and it shows that some aspects of its classical 
counterpart can still be preserved in the quantum vH model, and in particular 
the SG+FE phase. However, spin flipping introduced by the presence of T in 
the quantum vH model can modify the phase diagram, suppressing for instance 
the presence of SG+FE phase [2D]. However, it is well known that an additional 
transverse field in the KIL model with a Gaussian random bond coupling be- 
tween the localized Ising spins operators can produce important consequences 
as, for instance, a QCP [21] . 

In the present work, we will, therefore, study the KIL model with both the 
vH type of disorder for the intersite exchange interaction J y and a transverse 
field r which allows also to investigate the possible consequences for the phase 
diagram with the spin flipping. There is also another very important aspect re- 
lated to the vH type of disorder introduced in the present work: in the previous 
approaches using the Gaussian random bond SK-type J y [TpJ [TT] [12j HH [21] or 
the site disorder type given by the product of random variables the disor- 
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der is treated using the so called replica symmetry solution for the SG order 
parameters [52]. This solution is well known to have a serious flaw, because it 
is locally unstable below the freezing temperature [23- Certainly, that prob- 
lem could be overcomed by the use of replica symmetry breaking schemes |24j . 
However, this kind of scheme increases the number of order parameters in such 
a way that the search for order parameter solutions in the KIL model becomes 
extremely complicated. Nevertheless, that is not the only problem with the use 
of replicas to treat the disorder in the KIL model. There are also indications 
that the presence of one or other magnetic solutions could be dependent on 
the particular kind of replica symmetry breaking schemes |251 . By contrast, 
that is not the case for the disordered Jy given in the vH model (see following 
equation @). The disorder can be treated without the use of replica technique 
as demonstrated in the classical and quantum vH models [T9] [20] • Thus, the 
present use of the van Hemmen description of the disorder in the KIL model 
improves considerably the description of the Kondo-Spin glass- Ferromagnetism 
competition in disordered Kondo systems. It is important to remark that the 
present work is typically a mean field theory as in Refs. jTOJ [IT] [T2J [IT] [2JJ. In 
particular, the Static and saddle point approximations are used here. The use 
of the first approximation can be justified since our goal is mainly to describe 
phase boundaries as discussed in Ref. [21] ■ The saddle point method is in fact 
exact here, as a consequence of the long range nature of the vH coupling. 

This paper is structured as follows. In the next section, we introduce the 
model and calculate the corresponding thermodynamics. The following section 
is dedicated to discuss the numerical solutions of the saddle point equations for 
the order parameters and to derive the phase diagram. Finally, the last section 
is reserved to the conclusions. 



2 General Formulation 

The starting Hamiltonian in the KIL model is given by: 

H = 5>ft£ + eo5>( s + JK^p- dt + Sj t s+ 

ij,s i,s i 



, z (i) 
-y JijSfjSfj - 2r s fi . 

In Eq. 4JJ, S}. = |[nf t - nfj, S+ = §J t = (£+)t, Sf t = f^fa + f}^, 

§ di = 4t*^ S di = ( S ti) f > Ha =/i S fia, fit = d\ s d is where f} s (f is ) and d\ s (d is ) 
are fermionic creation (destruction) operators of / and d electrons, respectively. 
The spin projections are indicated by s or L 

The random coupling Jy in Eq. ([TJ is given as in the vH model by: 

Jij = Jf&Vj+Vitj) + jf (2) 

where £j and r\i in Eq. ([2]) are random variables which follow the bimodal 
distribution: 

P( x ) = h5( x -l) + 6(x + l)}. (3) 
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In Eq. S(x) is the Dirac delta function. As discussed in the previous section, 
the coupling Jij given in Eq. © is an infinite long range coupling which gives 
exact solutions in the thermodynamical limit for the saddle point approximation 
used below. 

The partition function is expressed within functional formalism using anti- 
commuting Grassmann variables (pi S (r) and ipi S (T) associated to the / and d 
electrons, respectively as [TU1 121) : 

Z = J DWili)D{<p*<p) exp [A VH +A K + A ] . (4) 

In the static approximation (SA) [T3] the actions in Eq Q are given as: 

u i,j (5) 
+^[(iu} + fid)Sij - Ptij\i£.(u) 

where, in the first term of Eq. ([5]), the chemical potential [if has been absorbed 
in e . 

Jk 



at « EE N-sH^-sH] EE 



(6) 



IS W JS 



with 



s h = \Yj£W*M' ( 8 ) 



The action AJ| a * is given in the mean field approximation (see Ref [10 .). In 
the remaining components of the action Aq and Ayjj, spinors are used : 

<pM=( ^51 ), *•(!»)=( til ) 0) 
and the Pauli matrices are given as usual by: 

; I) *-(; 7) *-(; -0- « 

We follow a procedure close to that introduced in Ref. [TU]. Therefore, 
the Kondo order parameter A CT ~ A = -^E] i' l ftjcr( u, ) ( Pi<r( UJ )) can be introduced 

in the partition function Z defined by Eqs. (UJ-©. Then, the ip fields are 
integrated and we obtain the following result: 

Z/Z% = exp(-2Np\\*)Z eff (11) 

where Z® is the partition function of free d electrons and 



z eff = f z?(w)ex P (^+EE£>)^>)t->)] 



(12) 
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with 



gjHu) = [(iw - pe )I + pTaz\5 i: 



f3 2 J 2 \ 2 



(iw + fJ,d)$ij - 0U 



-I. 



(13) 



In Eg. (|13p . we use the notation |A 2 | = A 2 and I_ means the unitary matrix. 

Introducing given by Eq. @, the action becomes composed of 

two terms: one randomic and the other one ferromagnetic. They can be rear- 
ranged to introduce SG and FE order parameters in Z e ft. The details of such 
calculations are shown in the Appendix. 

The free energy is, therefore, given by: 



(3F = 2/3J K X 2 - lim -^lnZ e// . 



(14) 



Using the saddle point solution for Z e ff (see Eqs. (([24 [) -(|25 [) )). the free 
energy in Eq. (fT4"]) becomes: 



/3F = 2/3 J K X 2 + 2pj qi q 2 + pj m 2 - ^ln(det G^^hj)). 



(15) 



The Green function G^j (w\hj) is given in Eqs. (l26|) -(l3~T |) . In order to proceed 
to the calculations, we use in the last term of Eq. (|T5|) the approximation 
introduced in Ref. [IU] which decouples the random magnetic field hj from the 
Kondo lattice. Thus, we obtain: 



lndet (G^HM) « ^E ln [dctr^HM] 



(16) 



with 



T^(u\hj) = [{iuj - /3e )I - ajij + /3Ta x ] 5^ 



3 2 J 2 X 2 



SkR u 



Y- 

N^(i(j + fi d )-/3e k - 



(17) 



Now, in the last term of Eq. (|15[) . we can use self-averaging property 
lEj/fe;Ci) = ((f(mO))e n - Therefore: 



45> 



where 



In 



(18) 



f*7 



<(/&»?)» 



to 



dtd V P(Z,ri)f(S,ri) 



(19) 



Then, by assuming that fid = and £o = 0, the free energy can be found as 
(iF = 2(3 J k \ 2 + 2(3J qi q 2 + f3J m 2 



1 



dx In 



0D 



cosh 



x + H 



(20) 



cosh^/-(a;- J ff) 2 +^ 2 J 2 A 2 
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Figure 1: Phase diagram T/J versus Jk/ J for Jo/ J = 1.6 and T/J = 0. 



with 

H = /3 V /[2J( W2 +^i) + 2J m] 2 +r 2 . (21) 

In Eq. (f2"0)) . the sums over the Matsubara frequencies and over k have 
been done in a way similar to Ref. [TU]. We have also used here the usual 
approximation of a constant density of states for the d electrons, p = for 
— D < e < D. The use of this density of states allows a direct comparison of 
phase diagrams obtained in this work with previous ones given in Refs. [101 HU 
[TB] . Finally, assuming that the probability distribution P(£, rj) = P(£,)P(rj), we 
can compute ((...)), in Eq. ([HI) using Eqs. © and (TlT)j) . 



3 Numerical results 

The coupled saddle point equations for qi, 52, m and A can be obtained directly 
from Eqs. (f2"D)) - ([2"Tj) . The numerical solutions for such order parameters allow us 
to obtain the following phases: (i) paramagnetism (PARA) given by gi = q>2 = 0, 
m = and A = 0; (ii) the SG phase given by q\ — qi ^ 0, m — and A = 0; 
(iii) the mixed phase (SG+FE) given by q\ = q% ^ 0, m 5^ and A = 0; (iv) 
ferromagnetism (FE) given by q\ = q2 = 0, m 5^ and A = 0; (v) Kondo state 
where only A is different from zero. For numerical results, D/J= 12 is used. 

Phase diagrams giving temperature T versus Jk (in units of J) can be built 
for several values of Jo /J and T/J. In Figure (fl]), such a phase diagram is 
displayed for Jo/ J — 1-6 and T/J = 0. For this case, in the large Jk region 
there is only one solution which corresponds to the Kondo state. When Jk 
decreases, the Kondo solution disappears. Actually, it is substituted by the 
magnetic solutions PARA, SG and FE which appear in that order when T is 
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Figure 2: Phase diagram T/J versus Jk/J for Jo/ J = 1.3 and T/J = 0. 



lowered. In Figure ([5]), we take Jo/ J = 1.3 and T = 0. This decrease of Jo /J 
from 1.6 to 1.3 does not affect the Kondo state, but changes a lot the magnetic 
solutions. In Figure ©, the solution FE is replaced by the mixed phase SG+FE, 
while the size of the region where the SG solution exists remains almost the same 
as in Figure (TT]). In Figure ©, the transverse field V is maintained equal to 
and we take an intermediate value Jo/ J = 1.4. As in the two previous cases, 
the Kondo state is not really affected in the large Jk region, but the region 
of the magnetic solutions in the phase diagram is again modified. Besides the 
existence of SG and SG+FE solutions, when the temperature is decreased, there 
is also an additional FE solution at much lower temperatures. In other words, 
in a small range of Jo/ J (1.3 < Jo/ J < 1.6), the phase diagrams present several 
scenarios concerning the existence of magnetic solutions. In contrast, the Kondo 
state is robust to such changes of Jo / J. 

Furthermore, a new situation is obtained when the transverse field is turned 
on, as can be seen in Figure (@J). For instance, for V = 1.0, the Kondo solution 
is obtained for a value of Jk / J a little larger than that found previously for 
r = 0.4 or r = and simultaneously, the range of Jk/J where the magnetic 
solutions are found is increased. Moreover, for such a decrease of T from 1 to 0, 
the transition temperatures between the magnetic phases are clearly depressed. 
But the most important feature observed with the increase of T concerns the 
magnetic solutions, because the SG+FE and FE phases disappear completely 
and it remains only the SG phase for a sufficiently large T value. 
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Figure 3: Phase diagram T/J versus Jk/J for Jo/ J = 1.4 and T/J = . 

4 Conclusions 

In the present work, the KIL model has been studied with assuming that the 
inter-site spin coupling J^- between localized spins is a random coupling given 
by the van Hemmen model as given in Eq. ([2])- It has also been added to the 
model a transverse field T which mimics a Heisenberg spin-flipping term. 

The results are shown in Figures (HJ)~ (SJ ■ For T = 0, they basically display 
two regimes when the strength Jk of the Kondo interaction is varied in units of 
the component J of the coupling (see Eq. @). In the first regime obtained 
for large Jk values, there is only the Kondo phase. In contrast, the second 
regime with only the magnetic solutions SG, SG+FE and FE exists when Jk is 
decreased. One important point is the order in which the magnetic phases are 
found when the temperature is decreased. For instance, the SG phase is found 
at higher temperature. Then, it can appear a SG+FE phase. The pure FE 
phase is found only at the lowest temperatures. It is also important to notice 
that the existence of the different solutions SG, SG -I- FE or FE depends on the 
strength of the ferromagnetic component Jo (given in units of J) of the coupling 
Jij, as can be seen in the Figures ([I])-©. When V is different from zero, the two 
regimes discussed previously are affected. While the Kondo solution needs larger 
values of Jk to be found, the magnetic solutions found at lower temperatures 
disappear rapidly when T is increased. 

It should be emphasized that the present approach using the Jy coupling 
given by the vH model yields two important improvements with respect to 
previous approaches. The first one concerns the use of the replica method which 
is not necessary here to generate the thermodynamics. This is an important 
improvement with respect to the previous approaches using the bond disorder 
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Figure 4: Phase diagrams T/J versus Jk/J for J /J= 1.4 and three values of 
T/J : 0, 0.4 and 1.0. The dashed, dotted and full lines are results for T/J = 1.0, 
T/J = 0.4 and T/J = 0.0, respectively. The critical lines for T/J = occur 
at higher temperatures than those ones for T/J = 0.4 and T/J = 1.0. In 
particular, for T/J — 1.0 there is no more SG+FE and FE solutions. 



given by the SK-like Gaussian random in the KIL model [ITJI I2TI [121 HHl IPT] 
or using the previous Mattis-like approach [TB]. For instance, the presence of 
magnetic solutions in these approaches is quite dependent on which particular 
scheme of replica solution is used, as explained in the discussion of Ref. [35] . 

As our present results suggest, the second improvement concerns the particu- 
lar kind of site disorder given by the vH model introduced in the KIL model with 
a certain range of Jq, which allows to obtain magnetic solutions SG, SG+FE 
and FE phases when the temperature is decreased. In that sense, the weak- 
ness of the approach proposed in Reference [16] is overcome and we are able to 
introduce here a mixed phase SG+FE. 

Thus, our present calculation using the van Hemmen site disorder can de- 
scribe Cerium disordered physical systems such as CeNii_ x Cu x or CePdi_ x Rh x 
alloys. In particular, Figures 3 and 4 can describe the phase diagram oiCeNii- x Cu, 
with Jk increasing with an increasing Nickel concentration, by explaining the 
Kondo behavior observed for x close to 1 and by proposing a good interpretation 
of the complicate magnetic behavior observed for smaller x values. These are 
indications that the use of the van Hemmen site disorder could be useful to de- 
scribe physical systems such as CeNi\_ x Cu x or CePd\- x Rh x alloys, although 
the low temperature phase is in these alloys a Kondo-cluster-glass followed by a 
disordered ferromagnetic one. However, it is important to notice that canonical 
spins have been used in the present work. This description is obviously not 
sufficient to capture the complexity of the cluster glass state. However, earlier 
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results for a mean field formulation of the cluster glass indicate that there are no 
essential differences between canonical spins and clusters of spins, as far as the 
phase boundaries are concerned 26 . One can, therefore, expect that most of 
the previous discussion concerning the sequence of magnetic orders as a function 
of Jk can be preserved even if the problem is formulated in terms of clusters of 
spins instead of canonical spins as it is done in the present work. 

On the other hand, we are presently working on a theoretical description of 
the Kondo-Cluster-Glass, by solving exactly the problem in a small cluster with 
n s atoms interacting between them by a disorder spin glass-like interaction. We 
have already solved the problem with only n s = 3 and a disorder intercluster 
bonding given by the Sherrington-Kirkpatrick interaction [27]. We think that 
the van Hcmmcn approach is easier to treat and we are presently working on 
clusters with a larger number n s , in order to have finally a more local description 
of the Kondo-Cluster-Glass observed in some disordered Kondo Cerium systems. 

In conclusion, we have to remark that our van Hemmen-Kondo description 
yields considerable improvements with respect to previous theoretical models in 
the two following points, the non consideration of the replica method and the 
problem of the mixed SG+FE phase. The validity of the van Hemmen model, 
which does not use the replica trick method, has been discussed in detail and 
it has been shown that this model is perfectly able to describe the spin glass 
experiments and that it is simpler than the other models for a mathematical 
treatment [55J [55]. On the other side, our van Hemmen-Kondo model gives 
with decreasing temperature a SG phase, a SG+FE one and finally a ferro- 
magnetic phase and the intermediate SG+FE phase is a real mixed phase with 
together non zero SG and FE order parameters. This model gives a good ac- 
count for the experimental phase diagrams of disordered Cerium systems, such 
as CeNii_ x Cu x alloys, and can be used to have a more local description of the 
Kondo-Cluster-Glass phase. 
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Appendix 

In this appendix, we present in details the procedure which allows to introduce 
the SG and FE order parameters in the problem. First, the random component 
of Jij given in Eq. ^) can be rewritten as: 



PJ_ 

N 




(22) 



~w 



N 



2/3J 



N 



£(^fx^f) 
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while the ferromagnetic one is 



N 



N 



E^ 



Jo 
N 



E^ 2 



(23) 



The last terms in Eqs (|2"2")l and ([2"3")l vanish in the thermodynamic limit. 

The Hubbard-Stratonovich transformation can be used to linearize the action 
A$jf. Thus Z eff in Eq. QJJ) becomes: 



Z eff = 



+ 00 



N 



2 f-\-oo p-\-oG p-\-ao 

dqi / dq 2 / dq, 

v J — oc J — 00 

2 



Nrh 



dmexv [ -~—{q\ + ql + qf) . } 



where the function A(qi, q 2 , q$, m) in Ea. ([2"4")) is: 



lnA(qi,q 2 ,q 3 ,m) 



A(q u q 2 ,q 3 ,m) = J D(M) exp I EEC^Vl^vM 



with: 



/? 2 J|A 2 



The random field in Eq. (|2"B1) is 



/ij = V2/^ + *£j92 + (Vj + 0)93) + y/2(3J m 

The saddle point solution of Eq. (|24|) gives: 



qi = i 



and 



'2/3 J ( qi +q 2 ) 
E v 5 ? > = V^J~om 



(24) 



(25) 

J. (26) 

(27) 

(28) 

(29) 
(30) 
(31) 



The symbol (...) is the thermodynamical average and i 2 = —1 in Eqs. (12"5|) - (13"TT) . 
The integral over the Grasmann fields can be performed in Ea. (|25[) . leading to: 



A( qi ,q 2 ,m) =exp ^]ln(det ^ (u\hj))^ 



(32) 
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